Identification of perturbation modes and controversies in ekpyrotic perturbations 
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If the linear perturbation theory is valid through the bounce, the surviving fluctuations from the 
ekpyrotic scenario (cyclic one as well) should have very blue spectra with suppressed amplitude for 
the scalar-type structure. We derive the same (and consistent) result using the curvature perturba- 
tion in the uniform-field (comoving) gauge and in the zero-shear gauge. Previously, Khoury et al. 
interpreted results from the latter gauge condition incorrectly and claimed the scale-invariant spec- 
trum, thus generating controversy in the literature. We also correct similar errors in the literature 
based on wrong mode identification and joining condition. No joining condition is needed for the 
derivation. 

PACS numbers: 04.20.Dw, 98.80-k, 98.80.Cq, 98.80.Hw 



I. INTRODUCTION 

The issue of scalar-type structure generated in the re- 
cently proposed ekpyrotic scenario is shrouded with con- 
troversies by two opposing camps [1~3] and [4-10]; for 
an introduction to the scenario, see [11]. The main point 
of [1] is that the dominating solution viewed in the zero- 
shear hypersurface (gauge) in the collapsing phase hap- 
pens to show a scale-invariant spectrum. However, this 
mode was identified in [5,6] as a transient mode in the 
subsequent expanding phase, thus uninteresting. In this 
work we wish to add some additional points to [6]. We 
will show that the same blue spectrum is generated even 
in the zero-shear gauge by identifying the mode relevant 
in the later expanding phase. Apparently, the same final 
observable spectrum should be derived independently of 
the gauge conditions used, and our result confirms it. We 
also point out that the possible scale-invariant spectra 
and others argued in [2,3,12] are errored by identifying 
wrong modes (often based on ad hoc joining conditions) 
which are transient in the expanding phase, thus irrele- 
vant. 

Before we embark on studying the evolution of struc- 
tures through bounce using the linear perturbation the- 
ory, we would like to state clearly the provisions we need. 
In [7] Lyth has clearly shown that the linear perturbation 
theory breaksdown inevitably as the model approaches the 
singularity in a singular bounce (if such a bounce is possi- 
ble at all, [13]), see also §VI of [8]. Somehow, this strong 
conclusion is unduely ignored by many authors [1,2]. If 
the bounce is singular we cannot rely on the linear per- 
turbation theory. Thus, in the ekpyrotic scenario and 
other bouncing models considered in this paper we will 
explicitly assume that the bounce occurs before the lin- 
ear theory breaks down, and this requires the bounce 
to be smooth and nonsingular. Although the authors 
of [1] claimed that the bounce in the ekpyrotic scenario 
to be singular, assuming a nonsingular bounce in such 
a scenario is legitimate, particularly if we consider the 



currently unknown physics near the bounce. Since the 
consequence of the singular bounce is clearly resolved in 
[7] {i.e., the linear theory fails!) investigating the remain- 
ing window with nonsingular bounce would be important 
to clearly resolve the remaining issue. 

In a single component fluid or field, the scalar-type 
perturbation is described by a second-order differential 
equation with two solutions (modes). In the large-scale 
limit (to be defined later) we can often derive a general 
asymptotic solution with two modes, see eq. (7). In an 
expanding phase we can identify clearly which ones are 
relatively growing (C-mode) and decaying (d-mode). If 
the initial condition is imposed at some early expand- 
ing epoch the decaying mode is transient in time, and 
naturally we are only interested in the relatively grow- 
ing mode. If wc introduce a collapsing phase before the 
early big-bang phase, however, the conventional growing 
and decaying classification can be often reversed. Still, if 
the large-scale conditions are met (and, of course, if the 
linear theory as well as the classical gravity are intact), 
the general solutions in eq. (7) remain valid through- 
out the transition. Thus, in our observational perspec- 
tive situated in expanding phase we are interested in the 
initial condition imposed on the C-mode, even if it was 
subdominating (relatively decaying) compared with the 
other mode when the initial condition was imposed. Al- 
though we made this point clear in [6] , in this work we will 
reinforce it by deriving concretely the C-mode initial con- 
ditions coming from the quantum vacuum fluctuations in 
the two gauge conditions used previously. In this way, we 
hope we could clear some of the controversies concerning 
ekpyrotic scenario and others in the literature. §11 and 
III are reviews. §IV contains our main results with con- 
sequences analysed in §V. We set c = 1 = ?i. 
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II. BASIC EQUATIONS AND GENERAL 
LARGE-SCALE SOLUTIONS 

Wc consider the scalar-type perturbation in a flat 
Friedmann world model supported by a minimally cou- 
pled scalar field. Our metric convention follows Bardeen's 
in [14] 



ds^ = -a^(l + 2a)d77^ - 2a^p^adr]dx" 



5S(l + 2<p) + 27,a|/3 



(1) 



and X = a(/3 + a7); an overdot and a prime indicate time 
derivatives based on t and 77, respectively, with dt = adrj. 
The background is described by 



jj2 _ 87TG 



r + v 



> + 3H^ + V^ = 0, (2) 



where H = -. The basic perturbation equations are 
[15,16]: 
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V = 0, v!' + 



(3) 

w = 0, (4) 



where z = a^/H, and 



V = a6(j)^, ips^ = ip- {H/4>)8(j) = -{H/4>)5(j}^, 

" = -'Pxl'^^ Vx = f- Hx- (5) 

ipstf) and ip^ are gauge-invariant combinations which 
are equivalent to the curvature perturbation (</?) in the 
uniform-field gauge {5(j) = 0, equivalently the comoving 
gauge) and in the zero-shear gauge (x = 0), respectively 
[14]. The perturbed action was derived in [16] 



(6) 



In the large-scale limit, meaning for negligible k"^ 
terms, eq. (4) has general solutions [15,16,8] 



ips<j,{k,r]) = C{k) - d{k) 
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To the higher-order in the large-scale expansion each of 
the four solutions have [1 -|- J2n=i 2 3 ^n{k\r]\)'^"'] factor 
with c„ differing for the four cases. We emphasize the 
general nature of these solutions in the large-scale limit. 
These are exact solutions of the spatial curvature pertur- 
bation (ip) in the respective hypersurfaces (gauges) valid 
as long as the fc^ terms in eq. (4) are negligible; thus 
valid for general (time- varying) potential V{(j)). Similar 
general solutions exist for the fluid situation for general 



(time-varying) equation of state P{lj), and even for the 
generalized gravity theories [17]. 

Schwarz has pointed out that as a smooth bounce has 
to violate the weak energy condition if space-time is flat, 
z"^ DC n + P will have (at least two) zeros; this means 
that 1/z will be ill defined and the higher order correc- 
tions will have singular coefficients, thus the long wave- 
length expansion becomes inconsistent, [19]. To achieve 
a bounce, in §V.C of [8] we have used an additional pres- 
ence of an exotic matter X with negative energy density. 
Thus, such an X-matter cannot dominate even during 
the bouncing phase. In [8] we have shown that if we con- 
centrate on the evolution of curvature perturbation, as- 
suming near adiabatic initial condition in the collapsing 
phase, eq. (4) for u, thus our solution for ip^ in eq. (7) as 
well, remains valid. In this context, z'^ goes through van- 
ishing points at least twice in the bouncing phase, and 
indeed, in that case the next order large-scale expansion 
includes / (l/z2)(i?7-order terms which are ill defined; one 
such term already appears in the d-modc of ips,j, in eq. 
(7) which is exactly the next order contribution. 

As parts of the series solutions in eq. (7) are ill de- 
fined for z = 0, in such a case wc should qo hack to our 
original equations (forms before we combine to make a 
second-order equation). One such original equation is 
conveniently available in the second equation of eq. (3) 
which shows that, for 2; = we have {(p-f^a/H)' = 0. 
Thus, for = we have an exact solution: ip^ cx H/a. 
Notice that our eq. (7) includes the above solution as a 
case! Therefore, we conclude that throughout the bounce 
(including z = points) , our leading order aymptotic so- 
lution for in eq. (7) remains valid. Thus, the ill 
defined higher order corrections in the series expansion 
do not cause any practical problem in the perturbations. 



III. POWER-LAW EXPANSION 

A field with an exponential potential supports power- 
law expansion/contraction of the scale factor [18] 
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In the power-law case eq. (4) leads to Bessel equations 
for V and u with different orders. Using the quantization 

based on the action formulation in eq. (6), we have the 
exact mode function solutions {p ^ 1) [20,6] 
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where x = k\ri\. The quantization condition implies 
|c2p — |cip = ±1 depending on the sign of r], [6]. 



IV. MODE IDENTIFICATION 

The Hankel functions can be expanded as [21] 

" -n^iy a) sinz^TT \2) T{v + n+l) 



n=0 
+ 
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T{-v + n+l) 



(10) 



Notice that, in the small x limit, the first (second) term 
in the parenthesis dominates for v < Q (y > Q). In eq. 
(7) the leading orders of the (7- modes are time indepen- 
dent whereas the leading orders of the d-modes behave 
as ips^ oc Iryp''" and ip-^ oc iTyp"^". Since 

^s^k cx \vr^Hll'^\k\ri\), cp^k cx k-'\T,\^-Hil'^\k\rj\), 

(11) 

wc can easily identify the first and the second terms in the 
parenthesis of eq. (10) as the d-mode and the C-mode, 
respectively. 

The lower-bounds of integrations in eq. (7) give rise to 
terms which can be absorbed to the other modes. Such 
an ambiguity is removed, for example, by identifying the 
C-mode of (fis4i in expanding phase in the large-scale limit 
when the time-dependent part of the d-mode has asymp- 
totically decayed away. The authors of [22] have shown 
that the d-mode does not necessarily decay away imme- 
diately after the horizon crossing; in some inflationary 
models wc have the d-modc effect not negligible near the 
horizon-crossing, and the final result can be interpreted 
as an amplification of the spectrum. Such an amplifica- 
tion occurs because, in expanding phase, it takes some 
time to have the time-dependent part of d-mode be negli- 
gible. The point is that there is no such an effect from the 
d-mode while in the asymptotically super-horizon scale. 
In our case of the bounce the relevant scale remains in the 
asymptotically super-horizon scale during the bounce, 
thus the solutions in eq. (7) arc well valid, and we do 
not anticipate any ambiguity arising while in the asymp- 
totically super-horizon scale. Accordingly, later in §IV, 
we will identify the C- and d-modes ignoring the contri- 
butions from lower-bounds of integrations of eq. (7). 

The power spectrum and the spectral index are defined 
as = ^\<fk\'^ and — 1 = dlnV^/dlnk. The spec- 
tral indices for the C-modes of (ps^, and ip^ can be read as 
(in the following, we assume the simplest vacuum state 
choice) 



{ns - l)v^^,c = {ns - l)v=^,c 



1-p 



(12) 



Although not interesting (because it becomes transient in 
an expanding phase) the spectral indices for the d-modes 
are 



{ns - l)^s<p,d 



6p 



1-p 



2p 
1-p 



■ (13) 



Notice that the spectral indices of the C-modc coincide 
in both gauge conditions, whereas the ones for the d- 
mode show strong gauge dependence. This is easily un- 
derstandable from the general solutions in eq. (7): in the 
power-law expansion case we have ipsif, /(p^ = 1 +p for the 
C-mode; the C-mode of (ps^, remains constant even under 
the changing potential whereas (p^ changes it value. Sim- 

ilarly, for the d-modc wc have Lpsci>/y^x ~ sp-i (^1^1)^' 
thus we have (ng - l)^^,d = {ns - l)^s4>,d - 4- 

We note again that the d-modcs show strong gauge de- 
pendence: the d-mode of ips4> shows more blue spectrum 
compared with (p^. Near singularity, the d-mode of (p^ 
diverges more strongly compared with the ones in the 
other gauge conditions [17]; see §111 of [8] for a summary. 
The strong divergence in the zero-shear gauge is known 
to be due to the strong curvature of the hypersurface 
(temporal gauge condition) [23]. According to Bardeen 
the behavior of ip^ "overstates the physical strength of 
the singularity", [23]. Thus, even in the collapsing phase 
where d-mode is the proper growing solution, we should 
not attach more meaning to the d-mode of <p-)^ than to 
the one of ips^. 



V. CONSEQUENCES 

In an ekpyrotic scenario with < p ^ 1 we have a 
very blue rtg — 1 ~ 2 spectrum for the C-mode. Although 
— 1 ~ for the d-mode of (/?^, we are not interested 
in the d-mode. Incidentally, we have ns — 1 ~ 4 for 
the d-mode of ps4> which better characterizes the physi- 
cal strength of the growing perturbation during the col- 
lapsing phase than Lp^. Our point is that, although the 
d-mode is the relatively growing solution in the collapsing 
phase, our classification of the C- and d-modcs is based 
on the general large-scale solutions in eq. (7). The large- 
scale conditions used to get these solutions are well met 
during the transition phase in the ekpyrotic scenario. In 
[8] we have shown analytically that, as long as the linear 
perturbation is valid, the solutions in eq. (7) remain 
valid throughout a (smooth and nonsingular) bounce, 
thus there occurs no mixing for the eventual growing so- 
lution in the later expanding phase. Thus, claiming the 
scale-invariant spectrum based on the d-mode of ip^ is 
incorrect; see the next paragraph for some technical de- 
tails. The C-modcs of both ips^ and p^ show the same 
blue spectra, and the complete spectrum of the C-mode 
and the one for the tensor-type perturbation can be found 
in [6]. 

We would like to comment on several minor compli- 
cations made in [1]. Firstly, the authors of [1] claimed 
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that by combining the scalc-invariancc of the d-modc of 
(p^ before the bounce and the new joining condition in- 
troduced by the authors they can derive a scale-invariant 
final spectrum. This implies that mixing occurs so that 
the rf-mode before the bounce sources and dominates the 
C-mode in the subsequent expanding phase. This con- 
tradicts with our result based on the general large-scale 
solutions in eq. (7). It was shown in [5,6] that the well 
known joining condition based on equations of motion 
[24] in fact confirms our result: i.e., in the large-scale 
limit the C-mode is affected only by the C-mode of the 
previous era. The new joining condition used in [1] is 
ad hoc and is not based on proper physical or mathe- 
matical arguments, see [6,9]. As we have shown in this 
paper, and more properly in [8], in order to trace the 
large-scale evolution of the eventual C-mode in expand- 
ing phase, wo can use the analytic solutions in eq. (7), 
thus we do not need the joining condition at all. Sec- 
ondly, in [1] it was emphasized that before the bounce 
the potential is restored to zero so that expansion rate 
changes to p ~ ^. As the perturbation still remains in 
superhorizon scale during the bounce such a change in 
the field potential does not affect the already generated 
perturbation spectrum. We have emphasized that the 
large-scale general solution in eq. (7) remains valid even 
under such a changing potential. We are interested only 
in the C-mode and the solution shows that the C-mode 
is not affected by the changing potential. Thirdly, the 
authors of [1] also stressed that radiation is present af- 
ter the bounce. In §V.C of [8] we have shown that the 
evolution of adiabatic (curvature) perturbation is not af- 
fected by the changing background equation of state or 
the presence of multiple component while in the super- 
horizon scale. Thus, the presence of radiation component 
after the bounce adds only a minor complication which 
does not affect the curvature perturbation in the super- 
horizon scale. 

In a similar context, for p = | we have ng — 1 = for 
the (i-mode of ips^f, (in this case we have ns — f = —4 for 
the d-mode of (f^). Identifying this as another possibility 
for generating a scale-invariant spectrum attempted in 
[12] is incorrect for the same reason as in the ekpyrotic 
case; this was pointed out in [6]. For the C-mode we have 
— 1 = 6, thus too blue. 

Another similar error was made in [3] , now in the case 
of p = ^. In this case we have ng — 1 = 4 for the C- 
modcs, and 2 for Lps,p and —2 for (p^ for the d-modes. 
Based on the zero-shear gauge authors of [3] claimed 
that the generated spectrum has = — 1 for P = 
and = 1 for p ~ (ekpyrotic!), both of which are the 
ones for the d-mode of (fi-)^, thus irrelevant for the final 
surviving (observationally relevant) spectrum. Although 
[3] used a bounce model which differs slightly from the 
one used in [8], as we have argued, while in the super- 
horizon scale the final surviving spectrum is independent 
of the changing background expansion rate. The authors 
of [3] also considered a radiation dominated era during 
the quantum generation stage whereas we considered a 



scalar field dominated era with p = i. It is well known 
that the scalar field with an exponential potential can 
be effectively identified as an ideal fluid with constant 
w{= P/l-i), thus the two systems coincide for p = ^. 

Yet another similar errors were recently added in the 
literature, [2]. The authors of [2] argued that one can- 
not ignore the entropy generation near the bounce of the 
ekpyrotic scenario; if the bounce is singular, we already 
have stated that the problem cannot be handled using 
the linear theory. Based on this argument the authors 
claimed that the conventional joining conditions should 
be changed. Unless we use the proper joining conditions 
derived in [24] , we can show that the growing (and dom- 
inating) d-mode in the collapsing phase can easily domi- 
nate and source the C-mode in the subsequent expanding 
phase while in the large-scale. In this way, the authors 
claimed ng = 1 spectrum for the ekpyrotic scenario which 
comes from the d-mode of (p^. However, we can see that 
the entropy generation anticipated near bounce would 
not affect the superhorizon evolution of (the C-mode) 
perturbation. The joining conditions known in the liter- 
ature give the same result as our present one based on 
analytic solutions, [5,6]. Perhaps the entropy generation 
would be important for the background evolution so that 
we could achieve a smooth and nonsingular bounce as we 
have investigated using toy models in [8]. 

The authors of [2] also have claimed, that even for 
the pre-big bang scenario the final spectrum should pick 
up the d-mode of p^ generated in the collapsing phase. 
For the pro-big bang scenario based on a conformally 
transformed Einstein frame wo have oq. (8) with p = i; 
thus we have a vanishing potential. In such a case we 
have — 1 = 3 for the C-mode, and 3 for ips^ and 
— 1 for (/3^ for the rf- modes. Based on the same logic as 
their ekpyrotic case, the authors claimed that the final 
spectrum should be rig = which is the one for d-mode of 
(fix- We already have explained what is wrong with such 
analysis and result. The correct = 4 spectrum in 
Einstein frame was derived in [25]. In the original frame 
based on the low-energy effective action of string theory 
the pre-big bang scenario shows a pole-like inflation with 
ace \t—to\~^^^. The perturbation spectrum was derived 
by us in the original frame with ns — 1 = 3, [26]. We 
also have shown that ipscf, is conformally invariant [27]. 
Thus, it is natural for the final spectra from the two 
frames (despite their very different descriptions of the 
background evolutions) to coincide. We note that in the 
orignal frame the pre-big bang scenario does not involve 
a contracting phase, and is just another inflation. In such 
a case, as in the case of ordinary inflation, the calculation 
does not require any joining condition. 

If the linear perturbation theory is valid throughout, 
and the bounce is smooth and non-singular (see [8] for 
several examples) we could rely on solutions in eq. (7) 
as long as the large-scale conditions are met. In such a 
case, as emphasized in [6], we do not need to use join- 
ing condition which, if we use the ones derived properly. 
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also gives the same result [5,6,8,9]. Several possibilities 
to have (smooth and non-singular) bounce models were 
studied in [8] . In [8] , using a toy bounce model based on 
an exotic matter with a negative energy density we have 
shown analytically that the pre- and post-bounce results 
of the C-modes of ^ipg^ and show the same behaviors as 
the ones we studied in this work (which ignores the pre- 
cise physics of the bounce) , independently of the presence 
of the exotic matter (and the boimce itself) introduced 
to connect the collapsing and the expanding phases, see 
§V.C in [8]. 

Equation (12) shows that the only way to get a 715 — 1 ~ 
spectrum from the power-law expansion based on an 
exponential potential is to have p » 1 which is the ordi- 
nary power-law expansion or a damped collapsing phase. 
As pointed out in [6], in the latter case as the model ap- 
proaches the bouncing phase the comoving scales shrink 
faster than the Hubble (dynamical) horizon. Thus, the 
large-scale condition can be violated near the bounce, 
and we cannot simply trace the perturbation through the 
bounce, see [8]. Therefore, the only remaining possibility 
to get an observationally viable spectrum is the former 
case which is just a well known version of inflation. 
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